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The zero-range potential approach is extended for the description of situations where two-body 
scattering is resonant in arbitrary partial waves. The formalism generalizes the Fermi pseudopo- 
tential which can be used only for s-wave broad resonances. In a given channel, the interaction is 
described either in terms of a contact condition on the wave function or with a family of pseudopo- 
tentials. We show that it is necessary to introduce a regularized scalar product for wave functions 
obtained in the zero-range potential formalism (except for the Fermi pseudopotential). This metrics 
shows that the geometry of these Hilbert spaces depends crucially on the interaction. 

PACS numbers: 03.65.Nk,03.75.Ss,05.30.Fk,34.50.-s 
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I. INTRODUCTION 
A. Zero-range approach in ultracold physics 

The roots of the zero-range potential approach he in 
the works on the neutron-proton interaction in nuclear 
physics. First, Bethe and Peierls JJ showed that it 
was possible to get numerous scattering properties on 
deuteron by replacing the neutron-proton force with a 
suitable boundary condition on the wave function in the 
s-wave channel. One year after, Fermi introduced 
the idea of a zero-range effective potential allowing a 
Born approximation in the computation of the scatter- 
ing cross section between a slow neutron and a proton. 
Equivalence between these two approaches was developed 
further by Breit Q and led to the final formulation of 
the so-called Fermi pseudopotential Q. The fact that 
the range of the nuclear force (~ 210~^^ m) is much less 
than the s-wave neutron-proton triplet scattering length 
(~ 4.310~^^ m) was a key ingredient for the success of 
this approach. 

More than half a century after these pioneer works, 
zero-range pseudopotentials appear as reference tools in 
the physics of ultracold atoms. There are several reasons 
for this renewed interest. First, the mean interparticle 
distance in ultracold gases is very large as compared to 
the range of interatomic forces (denoted by R), so that 
two-body collisions occur at energies very much smaller 
than the characteristic high-energy scale /mE? (m is 
the particle mass). Consequently, the full two-body po- 
tential contains useless informations (deep bound states, 
hard-core repulsion . . . ) for the study of the ultracold 
dilute gases. In this respect, the pseudopotentials which 
by construction reproduce the correct scattering phase 
shift at small energy, contain the minimal number of 
parameters for modeling two-body collisions. A second 
key property is that for processes of vanishly small en- 
ergies, a large class of pseudopotentials give the correct 
results at the first order Born approximation while this 
approximate treatment is not possible using the "real" 
two-body potential. Hence, mean-field schemes like the 
Gross-Pitaevskii and the Bogoliubov approaches |^ 



for the Bose gas or BCS 7j for the dilute two-component 
Fermi gas can be implemented using the Fermi pseudopo- 
tential. This idea has been generalized by introducing a 
family of pseudopotentials equivalent to the Fermi pseu- 
dopotential: the A-potentials 0. Varying the free pa- 
rameter A does not change the exact treatment of the 
two-body problem, while the Born approximation can be 
adjusted to the exact solution for a finite coUisional en- 
ergy. In the case of the ultracold Bose gas and using the 
Hartree-Fock-Bogoliubov formalism, the A-potential per- 
mits to self-consistently take into account the back action 
of the excitations on the condensate in three-dimensional 
systems Q and also in quasi-two-dimensional (2D) traps 
|9|. Scattering properties in low dimensional configura- 
tions are another powerful application of the zero-range 
potential approach. For example, in the case of s-wave 
scattering, the Bethe-Peierls boundary condition on the 
3D wave function is a very efficient way for finding the 
scattering amplitude of two particles trapped in planar 
or linear atomic waveguides _1Q, 11, 12J. The important 
concept of tunable interaction in quasi-one-dimensional 
and quasi-two-dimensional systems has been discovered 
thanks to this approach jlOl [ll| . The few-body problem 
is another area of research where the zero-range potential 
approach has already proven to be very fruitful. Exact 
results have been found for three and four particles in- 
teracting by pairs in the s-wave channel [isL IT^ ITsl [l6| . 
These studies are of first importance in the context of 
the resonant s-wave scattering regime obtained with Fes- 
hbach resonances these experiments the scatter- 

ing length can be tuned at will, so that it is possible to 
achieve the BCS-BEC crossover of the two-component 

Fermi gas El E mill 113. The zero-range approach 

permits us to evaluate the lifetime (which is remark- 
ably large in the neighborhood of the resonance) for the 
weakly bound molecules that appear in the BEC phase 
and also to find the value of the molecule-molecule scat- 
tering length in terms of the two-body scattering length 
|l4j . For the unitary quantum gas (infinite two-body 
scattering length) exact results have been found for the 
few- and many-body problems These studies 

show clearly that the zero-range potential approach is an 
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efficient tool beyond the two-body physics and permits 
us to obtain non trivial properties in the few- and many- 
body problems. 

More generally, study of strongly correlated systems 
obtained in the vicinity of scattering resonances is actu- 
ally one of the most challenging directions in the field of 
ultracold atoms. These systems are not only interesting 
in themselves but also permit accurate studies of issues 
raised in the quantum many-body problem. The BCS- 
BEC crossover in two-component Fermi gases is an exam- 
ple of such a situation 25, 26, 27] . The domain now en- 
larges to resonant scattering inp-wave (28, 29, SO. IsHIS^ 
and d-wave channels [s^ allowing possible studies of the 
BCS-BEC crossover in channels of high angular momen- 
tum in the near future. Major interest in these systems 
stems from the fact that the strength of correlations can 
be tuned arbitrarily while the mean interparticle distance 
remains large with respect to the range of interatomic 
forces R (for example, the s-wave scattering length can be 
adjusted to several orders of magnitude larger than R). 
Consequently, one can expect that short range physics 
is not directly involved in a large class of many-body 
properties which can then be described only in terms of 
the low energy two-body behavior. It is worth point- 
ing out that this fundamental feature has been verified 
in the BCS-BEC crossover of the two-component Fermi 
gas. For this reason, the zero-range approach which is 
parametrized only by the two-body low energy physics 
is a very appealing tool for studying these regimes. Ac- 
tually, the state of the art is as follows: broad s-wave 
resonances (like the one in ^Li at 835 G 20]) can be ac- 
curately parametrized using the Fermi pseudopotential, 
and narrow s-wave resonance js^] can be studied using 
a generalization of the Bethe-Peierls boundary condition 
|l5j ]. Recently, a general zero-range potential approach 
has been introduced for p-wave resonances and the 
aim of this paper is to extend this formalism for a de- 
scription of resonant scattering regime in arbitrary par- 
tial waves. 



B. Which generahzation of actual pseudopotentials 
is needed? 

The idea of replacing a true finite range two-body po- 
tential by a zero-range pseudopotential acting on each 
partial wave has been developed by Huang and Yang 
in the context of the hard-sphere model. However, due to 
the importance of the short range behavior in this system, 
pseudopotentials describing three-, four- and many-body 
correlations were also proposed for modeling the system 
beyond the dilute limit. The actual situation in ultracold 
atoms is radically different as the mean interparticle dis- 
tance is very much larger than the range of interatomic 
forces, and two-body pseudopotentials are sufficient to 
describe the low-energy behavior for the many-body sys- 
tem even in strongly correlated regimes. Moreover as the 
hard-sphere model cannot support resonances, the tools 



developed in Ref. 36] cannot be used for our purpose. 
Recently, a generalization of the Fermi pseudopotential 
for arbitrary partial waves and scattering phase shifts has 
been proposed using delta-shell potentials in the limit of 
small shell radius ,37.] , by the way correcting a mistake 
of Ref. Boundary conditions associated with the 

zero-range approach in Ref. |37| are defined for scattering 
states in free space and depend explicitly on their wave 
number k. Consequently, this approach applies only in 
time-independent situations. Instead in the present pa- 
per, we show that in resonant regimes, general (that is, 
energy independent) boundary conditions can be found. 
Another feature of the pseudopotentials in Ref. [s^ is 
that their action is defined on the radial wave function of 
the specified partial wave. However, it is of interest when 
dealing with wave functions of arbitrary symmetry (for 
example in the presence of an anisotropic external poten- 
tial) to have the full expression of the pseudopotential, 
that is, an expression which contains implicitly the pro- 
jection operator over the interacting channel (denoted in 
the following by H; for the channel of angular quantum 
number I). 



C. Organization of the paper 

The paper is organized as follows: in the first part, the 
regime of interest (two-body resonant scattering) is in- 
troduced and characterized in each partial wave channel 
by a set of two parameters in the scattering phase shift. 
In the second part, the different tools of the zero-range 
approach in this regime are derived. We show that the 
problem can be defined in terms of boundary conditions 
on the wave function which generalize the method intro- 
duced by Bethe and Peierls. We express these boundary 
conditions in spherical coordinates and also in Cartesian 
coordinates using the symmetric trace free tensors used 
in the usual multipolar expansion. This way of defining 
the zero-range approach permits us to include quite eas- 
ily within the expression of the pseudopotential both the 
projection operator H; and the Ith derivatives of the delta 
distribution. We show that in each interacting channel 
there exists a family of pseudopotentials generated by 
a free parameter, thus generalizing the A-potential ap- 
proach @,|3^. The third and last part, illustrates the 
formalism. We find the expression of the two-body scat- 
tering states from the pseudopotential and the Green's 
function method. By choosing a specific value of the pa- 
rameter A, this exact result is also obtained in the first 
order Born approximation for two-body processes at fi- 
nite energy. The pseudopotential of Ref. [33 is recovered 
as a particular limit of the present approach. Finally, 
the scattering states and also the low-energy bound states 
are found to be not mutually orthogonal in the zero-range 
scheme (except for the Fermi pseudopotential). We intro- 
duce then a regularized scalar product which solves this 
inconsistency. We consider the low-energy bound states 
supported by the pseudopotential and show that the new 
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scalar product leads to normalization constants which co- 
incide with the results given by the method based on the 
analyticy of the scattering amplitude . 



II. ZERO-RANGE FORMALISM IN THE 
RESONANT REGIME 

A. Regime of interest 

Let us consider two particles in the absence of an ex- 
ternal potential interacting via a short range isotropic 
potential U{f) (r being the relative coordinates between 
the two particles) . We denote their reduced mass by /x. In 
order to simplify the discussion and without loss of gen- 
erality, we work in the center-of-mass frame. The asymp- 
totic form of the scattering wave function ) with rel- 

ative wave vector k and coUisional energy E — 



2fi 



de- 



fines the partial amplitudes /; in each partial wave chan- 
nel of angular quantum number I through the relation: 



where kr ^ 1, ft = f/r, fik = k/k and Vi is the Legendre 
polynomial of degree I. The partial amplitudes can be 
also expressed in terms of the phase shift Si with: 



= -i- (exp(2*,5,) - 1) 
2ik 



(2) 



Typical situations in ultracold atoms corresponds to bi- 
nary collisions of particles having small enough relative 
velocities to ensure kR ^ 1, where B is the potential 
range. For two neutral alkali atoms in their ground state, 
the asymptotic form of U{r) is a van der Waals tail: 
U{r) ~ Cg/r^ and the criterion for low-energy scatter- 
ing processes (collisional energy E ^ En) is fixed by 

m 



1/4 



(3) 



Standard values for R (and En) range from 3.3 nm 
for ^Li {Er ~ 7.25 mK) to - 10 nm for ^^^Cs {Er ~ 
32.2 /iK). In this paper we consider situations where the 
phase shift can be parametrized in the following form: 



k 



2'+icot(5z(fc) = - — 



Wl 



aik 



(4) 



Eq.Q) corresponds to the so-called effective range ap- 
proximation for a short range potential |4l|: the param- 
eter Wl characterizes the Wigner threshold regime and 
—2ai is "the effective range" even if it has the dimension 
of a length only for I = 0. For neutral particles interact- 
ing with a potential of asymptotic form U{r) ~ Cn/r^ 



(r oo) the threshold behavior is valid for I < (n — 3)/2 
and the effective range approximation for I < (n — 5)/2. 
In ultracold atoms in their ground state the power law is 
given by n = 6 so that the expansion in Eq. is rigorous 
only in the s-wave channel (wq is the scattering length), 
while in the p-wave channel only the first term in the ex- 
pansion is justified. For I > 2, even the Wigner threshold 
regime given in Eq.l^J is not valid. However, we adopt 
this form as a generic way to parametrize the resonant 
regime in an arbitrary partial wave channel. 

In broad s-wave resonances, the effective range is neg- 
ligible (ao — 0) and \wo\ is very large as compared to 
the potential range R. For wq > 0, /o has a pole as- 
sociated to a shallow bound state of energy —h^ /2fj,WQ. 
Its wave function coincides (for r > R) with the molec- 
ular state populated by pairs of particles in the BEC 
region of the BEC-BCS crossover. Remarkably, in the 
neighborhood of the resonance, the scattering cross sec- 
tion takes large values whatever the sig n of the scatter- 
ing length Wq. As explained in Ref. |38l |. in the channels 
I > 0, this property is no longer verified, moreover the 
"effective range" parameter is essential for a modeling of 
the shape of the scattering amplitude and of the shal- 
low bound state both. To define the resonant regime for 
Z > 0, we follow the discussion given in Ref. and use 
the expression of the partial amplitude associated with 
the phase shift in Eq.Q): 



Mk) 



wik"^ 



1 + wiaik^ + iwik^^~^^ 



(5) 



The resonant regime corresponds to situations where 
the parameter \wi\ takes arbitrary large values, so that 
\wiai\ ^ R^ . It can be shown [Sg that describing the 
resonant regime by Eq.© is valid only for a/ > (this 
assumption emerges also as a consequence of the regu- 
larized scalar product that we introduce in the last sec- 
tion). The main difference as compared to s-wave broad 
resonance is that the resonant character of the scatter- 
ing cross section associated with Eq. Q depends crucially 
on the sign of wi. This property can be easily seen by 
writting Eq.(|SJ) in the Breit- Wigner form: 



fi 



1 



r/2 



kE - Er + ir/2 



(6) 



with Er = 



— and r = . For large and 

2fiaiwi fiai 

positive values of wi, Er < and the scattering ampli- 
tude has a pole which corresponds to the presence of a 
real bound state of vanishing energy Eb ~ Er- In this 
regime, does not present a resonant structure. How- 
ever, for large and negative value of wi, the bound state 
transforms into a long-lived quasi bound state which pro- 
duces a resonance in the scattering cross section with 
~ l/\wiai\ in the neighborhood oi E — Er. This 
feature is a consequence of the existence of the centrifu- 
gal barrier in the radial equation for I > (see Fignj. It 
is worth pointing out that the resonant regime occurs for 
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FIG. 1: Schematic representation of a two-body effective ra- 
dial potential Vi(r) = [/(r) H J' ~!" — ^ in a partial wave 

i > 0. Er is the energy of a quasi-bound state 0, E3| which 
leads to a resonant behavior in the two-body scattering for a 
coUisional energy i5 ~ _Er. 



relative wave vectors k verifying kR <tc 1, which justifies 
a zero-range potential approach. 

In actual experiments on ultracold atoms, the reso- 
nant character in the two-body scattering is ensured us- 
ing Feshbach resonances ■ Interaction between atoms 
depends on the spin configuration of the two-body sys- 
tem and Feshbach resonances appear as a result of the 
coupling between a closed channel (supporting a short- 
range molecular state) and an open channel (associated 
with the scattering states). The energy of the short-range 
two-body molecular state can be tuned by use of an ex- 
ternal magnetic field. The scattering resonance occurs 
when the energy of the molecular state is in the vicin- 
ity of the collisional energy of the two scattering atoms. 
The molecular state in the closed channel depends on 
the short range properties of the system and is not de- 
scribed by the zero-range approach. Consequently, the 
bound state or quasi-bound state which appears in the 
zero-range model is associated only with the tail of the 
molecular state resulting from the coupling between the 
closed and open channels. 

In the s-wave narrow Feshbach resonance (large and 
negative effective range) of ^'^Na at 907 G [SJI, ao — 
26 nm, and R ~ 4.5 nm. In higher partial wave chan- 
nels, scattering is not isotropic i.e. it depends on the 
angular quantum number m of the state. For example 
the p-wave resonance in ""^K (2^] (i? ~ 6.8 nm) has been 
parametrized with the law in the states \m\ = 0, 1 as 
a function of the external magnetic field. There are two 
zero energy resonances (Er — 0) very close to one each 
other at 198.84 G (for m = 0) and 198.31 G (for \m\ = 1) 
with a width characterized respectively by aiR ~ 2.74 
and 2.81. Generalization of the zero range approach for 
taking into account this anisotropy adds technical de- 
tails not essential in a first approach and in this paper, 



the discussion is centered on the modeling of isotropic 
interactions. 



B. Defining tiie zero-range approach 

In this section, the interacting potential U{r) which 
can be neglected in the Schrodinger equation for |r | > i? 
is replaced by a specific behavior of the wave function 
as the relative coordinates between the two interacting 
particles goes formally to zero. This behavior defines 
boundary conditions on the wave function which depend 
on the parametrization of the scattering phase shift in 
Eq.gj). 

Outside the potential range (r > R), the scattering 
wave functions {'^^{r)} of the two-body problem are so- 
lutions of the free Schrodinger equation with the asymp- 
totic behavior given by Eas. (|ll2l4() : 

oo 

'^kir) = Y.T^i{n-nk)Rki{r) , (7) 

/=o 

where Rki (r) are the radial functions which can be ex- 
pressed in terms of the two spherical Bcsscl functions 
as ,43,j : 

Rki (r) = JVki ( ji (kr) - yi (kr) tan 6i{k)) . (8) 

The zero-range approach amounts to considering the for- 
mal continuation of Eq. Q for < r < i? . As a con- 
sequence, using the behavior of ji{z) and yi{z) at the 
vicinity of the origin: 

(2^-1)!! / zV2 A , ^ 

one can deduce the scattering states from the following 
boundary conditions of the wave function at r = 0: 

r™ 1^ (2/)!! tan<5z(fc)y' "'^ ^ 

This expression is a first step toward the formulation of 
the zero range approach, but does not have a general 
character in the sense that the quantum number k ap- 
pears explicitly. Instead, we search for a boundary con- 
dition which applies on any linear combination of scat- 
tering states with different collisional energies. This way, 
the formalism can be used when the collisional energy 
is not defined, as it can happen in time dependent situ- 
ations. Using the parametrization of the phase shift in 
Eq.Q), one obtains from Ea. ltTT)) : 

lim ( 9f + + 1 V'+ii?,,(.) = 0. 
r^o \ (2^)!! wi J 

(12) 
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The idea is then to substitute the energy dependent term 
in Ea. (|12|l by the appropriate radial derivative. From 
Eas. Hinil2|l we find the desired condition: 



wi 



where we have introduced the differential operator: 



{2l-l)aidt 



(13) 



(14) 



We are now ready to define the zero-range approach in 
general situations. The wave function ^E* of the interact- 
ing two-body system is expanded in each partial wave 
channel in terms of the spherical harmonics as: 



(r|n,|^') = 



I 

E 

m— — l 



(15) 



where H; is the projector on the subspace of angular 
quantum number and (r, fi) are the spherical coordi- 
nates with f2 = (6*, (/>)■ For r 7^ 0, \l/(r) is a solution of 
the free Schrodinger equation (it implies that the multi- 
polar radial functions c^{r) are regular functions of the 
interparticle distance r), and for an isotropic interaction, 
the potential term is replaced by the following boundary 
conditions on the wave function at r = 0: 



lim \Vi + — 

^0 V wi 



4""{r) = 



I < m < 



(16) 



Eq. (|16|l generalizes the Bethe-Peierls approach which 
applies in the s-wave channel when the effective range 
can be neglected (ao = 0): 



5,(r*) 
lim — — 



1 

Wo 



(17) 



where the limit \wq \ S> -R is associated with a zero energy 
broad resonance. 

Note that the generalization of Eas. ljltill to the case 
of an anisotropic two-body interaction is straightforward 
by introducing parameters that depend on the quantum 
number to: {'Wijn,aim}- 



Cartesian representation of the multipolar 
radial functions 



While the zero-range approach is now established in 
a general way by Eas. (|14ll5ll6() . it can be easier to use 
equivalent boundary conditions in the Cartesian coordi- 
nate system (for example when the system is inhomoge- 
neous or anisotropic). Furthermore, as will be shown in 
the next section, using these conditions leads to a simple 
construction of the zero-range pseudopotential (in par- 
ticular, derivatives of the delta distribution are easily 
introduced in Cartesian coordinates). In the following. 



I 









1 


Or — aodr 


1 


n° 


\dl + a^dl 


2 


n"n^ - i^"" 




3 


n°'nl^n< - \ {n°'5l'^'< + permut.) 





TABLE L Examples for the symmetric traceless tensors Si 
Eg. 1181 and for the differential operators Vi in Eg. 1131 . n" 
are the components of the unit vector n = r/r. 



we adopt the standard assumption of implicit summa- 
tion over repeated indexes, the three directions of space 
{x,y,z} are labeled by a Greek index a £ {1,2,3} and 
[ck/S . . . ] denotes a set of I indexes. The idea is to use the 
symmetric tensors Si which appear in the multipolar ex- 
pansion used in electrostatics p^. 



They can be defined 



through their components by: 



(18) 



where we have introduced a short-hand notation for the 
l-ih partial derivatives: 



d 



a/3.. 



_d d_ 



(19) 



and {2;"} are the three components of the radial vector 
r in the Cartesian basis {ca}: i.e. r — x'^Ba. These ten- 
sors can be expressed in terms of the component of 
the unit vector n = f/r and are functions of the angles 
Q, only. The normalization factor in Ea. ll8|l has been 
chosen for later convenience and explicit expressions of 
Si and 2?; with / < 3 are given in Tab. nil. In the follow- 
ing, we will use some useful properties of these tensors 
which are collected in the Appendix. By construction. 

Si are traceless {5af}S\"'^"'^ = 0) and symmetric tensors 
of rank /. They are eigenvectors of {L is the angular 
momentum operator) with eigenvalue Ti'^l{l + 1). Simi- 
larly to the spherical harmonics F'™(ri), the tensors Si 
can serve as a basis for an alternative multipolar expan- 
sion of a wave function, like in Ea. (|15|l . To this end, we 
define the Cartesian representation of the multipolar ra- 
dial functions of order I associated with a wave function 
^, by a symmetric trace free tensor C*_i(r) having the 
components: 



C, 



Q/3. 



An 



(20) 



For example, the monopolar (/ = 0) function C<s,fl{r) 
is just r times the projection of the wave function in 
the s-wave channel. Quite generally, C^f"'\r) are regu- 
lar functions of r analogous to the coefhcients c^{r) in 
Ea. (|15|) and characterize the behavior of the wave func- 
tion ^'(r, fl) in the channel I through the expansion (see 
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Eas. H69l72() in the Appendix) 
(2/ + 1)!! 



(rin^i*) 



(21) 



Using these notations, the generahzed Bethe-Peierls 
boundary conditions Ea. l|16|) can be written as: 



lim(A + -)4t-l(0=0 

r^O Wl 



(22) 



Eas. H2ni21l22(l are equivalent to Eas. p5ll6|l . We show in 
the next section that they permit to find a general form 
of the pseudopotential. 

Note that an alternative formulation of the condition 
(|22|l in the s-wave channel amounts to imposing the fol- 
lowing singularity on the two-body wave function: 

(f inolV-) = (l - - + Ao(r2 + 2aor)) + 0{^). 

r \ a / 

(23) 

In the p-wave channel I — 1, the short-range behavior 
of the interacting two-body wave function takes also a 
simple form: 

{f\Uim = ^ (l - ^ + - 2air3)^ + 0{r'), 

(24) 

where the "dipolar momentum" p defined by p" — 
3C§ i(0), characterizes the p-wave singularity of the 
wave function. In general situations the coefficients 
{Ao,Ai,Cqjfi{0), p) can be functions of time and of the 
center of mass coordinates of the two particles. 



D. Zero-range pseudopotential 

Instead of giving a final expression of the pseudopo- 
tential directly, we prefer to detail the different steps in 
its derivation. Such a way of presentation can be useful 
if one wants to derive pseudopotentials associated with 
others expressions of the phase shift in Eq.Q) or in dif- 
ferent contexts (low dimensions for example). In either 
cases, a zero-range pseudopotential can be deduced from 
the three following constraints: (i) it cancels the delta 
terms arising from action of the Laplacian on the singu- 
larities of the wave function at r = in the Schrodinger 
equation, {ii) it imposes the correct boundary conditions 
and (iii) it can be used in the Born approximation. The 
first condition (item-i) can be fulfilled by introducing po- 
tential terms in the different channels I of the form (see 
Eq.(|77I) in the Appendix): 



(-1) 



c 



\^){dU...]5W) ■ (25) 

In order to include the boundary conditions Eas. (|22|) 
(itcm-ii), we consider the "boundary operator" Bi de- 
fined by: 

(r 1^,1*) = (a^..]<5)(r ) lhn(P, + i-)4"_f-l(r). (26) 



By construction wave functions belonging to the kernel 
of Bi satisfy the desired boundary conditions. Conse- 
quently, the family of pseudopotentials defined by: 



V, 



f3Bi 



(27) 



where f3 is an arbitrary parameter, can be used to solve 
the two-body problem in the resonant regime. Interest- 
ingly, while exact zero-range solutions can be obtained 
V/5 ^ 0, it is not hard to be convinced that a first or- 
der Born approximation performed for example in the 
computation of the two-body scattering states leads to 
a /3-dependent result. Indeed in this approximation, the 
two-body wave function ^''^"^ is an eigenstate of the non- 
interacting Schrodinger equation (a plane wave in the ex- 
ample chosen) and is then a regular wave function of the 

radius r. Consequently, the action of Vj^*'' on gives 

a zero result while (r\vl'^ + PBil'^^^^) is /3-dependent. 
Furthermore, and we will show this property explicitly 
in the next section, varying f3 offers the freedom to ad- 
just the pseudopotential in such a way that the first order 
Born approximation is exact at a given coUisional energy 
(itcm-wi). Finally, the expression of the pseudopotential 
in Eq. 127|l although general, can be rewritten in the more 
convenient following form: 



(28) 



In this expression, g\ i = (— 1)'/3 is the coupling constant 
defined by: 



9x,i 



2TTi2l + l)h'^wi 
//(I - Xwi) 



(29) 



and TZx^i [ . ] are symmetric trace free tensors which gener- 
alize the regularizing operator introduced in Refs.,&..35j|. 
They act on a wave function 5* as: 



7^ 



a/3. 
A,/ 



(30) 



The parameter A appears in Ea. l29|) and Ea. H3U|l both, 
and by construction (like /? in Ea. H27|l ) it is a free param- 
eter of the theory: exact wave functions obtained in the 
zero-range formalism do not depend on the value of A. 
The pseudopotentials in Ea. (|28|) generalize for arbitrary 
partial wave channels the s-wave and p-wave A-potentials 
already introduced in Ref . ^ i35|| . 



III. SOME ILLUSTRATIONS AND 
PROPERTIES OF THE PSEUDOPOTENTIAL 



Two-body scattering states from the Green's 
method 



This section is both an illustration of the formalism 
and a way to test the proposed form of the pseudopo- 
tential given in Eq. (|28|) . We consider the case where two 
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particles interact in the Zth partial wave channel with- 
out external potential and we derive the scattering states 
from the general Green's function method. The interact- 
ing two-body wave function ^i:{r) (wave vector k, energy 
E = ?i^fc^/2/i) is then the solution of the integral equa- 
tion: 

^f^ir) = ^f{r)-jd^r'GE{r,r'W' \VxA^^) , (31) 

where ^^~\'r) = exp(ifc.r) is the incident plane wave 
and Ge is the one-body outgoing free Green's function 
at energy E: 

GE{r,r')^ expzfcz. ^.^^^ u=|f-f'|. (32) 
Integrating Eq.JSIJ over r ' coordinates gives: 

where 9'^^ j stands for the /th partial derivatives with 
respect to r ' coordinates. Thanks to the trace free prop- 
erty {5apT^xi ' \'^k\ ~ contribution of the Zth deriva- 
tive in Ea. H33|l takes a simple form. Indeed, for a given 
function F{u), one can show that: 

{n^np . . . )<f-l[*£](-r)' (^h^ F{r) , (34) 
and with F{u) = G£;(r, r'), Ea. (|^ reads finally: 



/ exp(iA:r) 



(35) 



As r — > 0, one recognizes a multipolar singularity in the 
two-body wave function: 

_ (2/ + l)!!^;5;;„g...] t,^...] 
^fein (i_Au;0 r'+i ^-^ + ^' 

(36) 

and comparison with Ea. ()21|l gives the multipole compo- 
nents: 



4tr'(o) 



Wl 



1 - Xwi 



n 



a/3.. 



(37) 



In the partial wave channels of angular quantum number 
different from I, components of the wave function '^i:(f) 

coincide as expected with the partial waves of *?^(f). 
All the information on the interaction are gathered in 
the term 5^[c</3...]7^^"f■■■'[*fe] in Eq.JSS). From Eq.^, 



one can see that it contains implicitly the projector 11; 
through: 



1 



limiVi + Xy+\r\Ui\^) . (38) 



{2l + l)\\ r^o' 

A closed equation is obtained by applying the operator 
5^[a/3...]7^[^"f ■■■'[■] on both sides of Eq.ljSHl. Using the 
usual decomposition of a plane wave over spherical ones: 

oo 

exp(ifc.f) = J2i'{2l + l)Vi{n.nk)ji{kr) , (39) 
together with Eq.Q this leads to: 

i2l + iy.\{ikyViin.nk) , (40) 



lim(A +A)r'+i (r\ni\^^y>) 

r—>0 fc 



on the other hand: 



lim(r>/ + A) 

r— >0 



,2/+l 



-A 



f exp(ifcr) 



r / \ r 



(-1)'(2Z-1)!! (A + aik' + ik''+') , (41) 
and finally, one gets: 



<? -pM-lr.T, 1 - ~ \wi){ikyVi{n.nk) 



(42) 



As expected, the expression of the scattering wave func- 
tion does not depend on the parameter A: 

''^ ' k ' 1 + wiaik^ + iwik'^i+^ 

(-.fc.)'(ia.)'(^) . (43) 

Asymptotically for kr ^ 1, the two-body wave function 
is given by: 



*^:(f ) ^ exp(ifc.f ) + {2l + l)riin.nk)fi 



exp(z/cr) 



, (44) 



where the partial amplitude /; is given by Eq.©. 

It is interesting to note that the exact result can be also 
obtained at the level of the first order Born approxima- 
tion (which is a "A-dependent" treatment) by choosing a 
particular value of the parameter A: 



2/-I-1 



(45) 



For example by setting A = in a first order Born approx- 
imation gives a correct evaluation of a scattering process 
at a vanishly small collisional energy (Wigner threshold). 

Generalization to the case where the interaction oc- 
curs in several partial wave channels is straightforward: 
it is obtained after simple addition of the corresponding 
pseudopotentials given by Ea. H28|l . 
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B. Links with other formalisms 

It is instructive to compute the matrix element of V\^i 
between two states |$) and IvI/). A simplifying step is to 
realize that the projection operator in the partial wave 
(I) is implicitly present in the expression of the pseudopo- 
tential Ea. H28|l . Indeed, assuming that $(r ) is regular at 
r = 0, one obtains: 



:/3...]*)-o = Z! lim 



(f|n,|$) 



(46) 



Moreover from Eas. H38l4(j|l and Ea. (|73|) in the Appendix, 
one finds: 

lim(A + A)r'+i(f|n,|*) ] . (47) 

r— >0 / 



First, we consider plane waves {f\k) = exp(ifc.r). The 
matrix element between two such states is: 



{k'\Vxj\k) ^ gx.iikk'Y Viirik' .Hk) 



(48) 



with fik' — k'/k'. For wave vectors of same modulus 
{k = k' = ^/2f^E/n) and for the particular value of the 
free parameter A given in Eq. H45|) , comparison with the 
expression of the partial wave amplitude Eq.(|Sll shows 
that Ea. H48|) coincides with the on-shell T-matrix t{k ' <— 
k) taken at energy E: 

t{k'^ k) = (fc'|F^Bo,.„,|fc) = -M!/(fc'^ fc), (49) 

where f{k' ^ k) = {21 + l)Vi{nk' .nk)fi{k) is the full 
scattering amplitude. 

Another interesting application of Eq. H47|) is given by 
considering spherical outgoing waves: 



vector kp- However, in this mean field approach the ap- 
proximation made on the interaction term is only valid 
for binary processes of vanishly small colliding energies 
(equivalent to a Born approximation with the pseudopo- 
tential and the particular choice A = 0). Consequently 
in the resonant regime where fci^jwol ^ 1, the character- 
istic coUisional energy of pairs of particles on the Fermi 
surface is /itiwq and the approximation made is not 
justified. By the way, it is possible to build a qualitative 
and simple mean field model which supports the stability 
of the gas at the BCS-BEC crossover 46] . 

Taking two states with angular parts characterized by 
the same spherical harmonic in Eq.|07j), one obtains the 
expression of the pseudopotential in radial coordinates: 



lim(P; + A)(r 



(52) 



In the particular case A = 0, a; = 0, r.h.s. of Ea. H52() 
coincides with the form of the pseudopotential introduced 
in Ref. El: 



».^5^^?iii^^liin3-'(r'«0 ^ (53) 



{r\klm)^ Y^"'{n)ji{kr) 



(50) 



2^(2/)! 



Note that the delta-shell regularization procedure intro- 
duced in Ref. jH3 while correct appears not to be essential 
in the zero-range approach. 



C. Regularized scalar product 

Except for the case of the Fermi pseudopotential which 
is a particular case of Ea. l|^ . two scattering states '5^, 

and '^^ (with k ^ k') of expressions given by Ea. H43(l 
are not orthogonal each others with respect to the usual 
scalar product. To examine this point, we exclude the 
singularity at r = in the scalar product between the 
two different scattering states by introducing a cut-off tq. 
We obtain from the Schrodinger equation: 



In this case, the matrix elements are given by: 



{klm\Vx.i\k'l'm')=5iv5-n 



gx.ijkk 
'47r(2/-M) 



(51) 



For A = 0, this result coincides with the one obtained in 
Ref. Here, the A-frcedom is a new ingredient and 

shows clearly that mean field calculations in zero-range 
approaches have to be developed carefully. For example, 
considering a two-component Fermi gas close to the uni- 
tary limit with — Q and a large and negative value of 
the scattering length wq, the mean field model in Ref. 
predicts an instability which is not observed in experi- 
ments mmmiHllll. TWs failure can be explained 
as follows. The interaction concerns only particles close 
to the Fermi surface (Pauli blocking) with a Fermi wave 



d3f^'l,(f)*^-(f) = 



r>ro 



d'^nmXr)dr-^dr) -■^dr)dr%,{r) ■ (54) 



This relation is general i.e. it is also valid for scattering 
wave functions associated with a finite range two-body 
potential whatever the value of rp . In standard situations 
where the Hamiltonian is hcrmitian, the r.h.s. of 1)54(1 is 
zero for tq = 0. However in the zero-range approach, 
for wave functions satisfying the boundary conditions in 
Eq. (|13() , as ro ^ this term converges toward a constant 
value for I = and ao 7^ and even worse, it diverges for 
I > 0. Interestingly, the regular part of Ea. H54|) for such 
wave functions can be obtained through the following 
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surface integration: 



Particles interacting in the s-wave channel only: 



Reg / d'f^l,ir)'^j:{f) 

oo 21+2 



d'nm,{r){r\Ui\-9t)A55) 



where the operator Reg { . } extracts the regular part at 

To = 0. The meaning of this result is that the Lapla- 
cian operator is not Hermitian with respect to the usual 
scalar product when one considers singular functions like 
the ones in Ea. H43|l . In another hand, the true scatter- 
ing states associated with the finite range potential ex- 
perienced by the particles, are orthogonal to each other 
so that direct identification of the true scattering states 
with the states {l^fj:}} is not possible. This feature fol- 
lows from the fact that the mapping between the two 
eigen basis is not valid for r R while it is justified only 
outside the potential range. Indeed, the singular bound- 
ary behavior imposed on wave functions in Ea. l|lt)|) is a 
way to reproduce the effect of the true finite range poten- 
tial for r > R but has a formal character for interparticle 
distances r ^ R. In order to understand this issue more 
closely, we consider the sphere of radius R which delim- 
its the outer parts (region r > R) from the inner parts 
(region r < i?) of the true wave functions. Contribution 
to the scalar product coming from the outer-part of the 
true wave functions is given by Ea. (|54|l with r^ — R and 
due to orthogonality between the true scattering states, 
scalar product of the inner-parts cancels this term. How- 
ever in the zero-range potential approach the inner parts 
of the wave functions (r < R) are not properly described, 
so that compensation of the two contributions does not 
occur. To solve this inconsistency, the idea is to modify 
the usual scalar product itself by including implicitly the 
contribution of the inner parts. We define then a regu- 
larized scalar product ( . | . )o by subtracting the surface 
term {r.h.s. of Ea. H54|) ') from the usual scalar product as 
ro^O: 



(* |$)o = Reg 



<iV**(f)$(f) 

r>ro 



oo 



21+2 



\.i™o^((2Z-l)!!)V.= 



d'n'f*{f){f\Ui\<i>), (56) 

1"0 



so that by construction, two eigenstates of the zero-range 
Hamiltonian are orthogonal to each other with respect to 
this metrics. 

In general there are several possible representations of 
Eq. (|56|l due to the various and equivalent ways of extract- 
ing a regular part. We give below explicit expressions of 
the regularized scalar product associated with the pseu- 
dopotential (|28|l in the I = and I = 1 partial wave 
channels: 



(*|$)o = lim i / dV^'*(r )$(f ) 

+aorl [ d^n^*{f)^r) } . (57) 



This result can be written more formally by intro- 
ducing a scalar product of the form: 



(*|$)o = ydV5(r)^'*(r )$(f) , (58) 

with the weight g{r) = 1 + aoS{r) . 

• Particles interacting in the p-wave channel only (for 
example in the case of polarized fcrmions): 

(^'|$)o = lim j / dV**(f)$(f) 

+ (air*-r3) / d^n^*{r)^r)\ . (59) 

This scalar product has been first introduced in 
Ref. jH3 in the same form as in Ea. l58|l with the 
weight g{r) = 1 + 5{r) [(air^ — r) . 

D. Bound states and their normalizations 

In this section, we consider the normalization of the 
two-body bound states of vanishing energy supported by 
the zero-range pseudopotential. We show that using the 
regularized scalar product is equivalent to the method 
based on the analyticy of the scattering amplitude [38l| . 
This last method was used in Ref. [T^ in the context 
of s-wave narrow resonances. In this case or also for p- 
wave or higher partial waves channel resonance where 
bound states are even not square integrable with respect 
to the usual scalar product (they don't belong to the 
Hilbert space), the regularized scalar product appears 
as the natural way to perform the normalization in the 
configuration space. Moreover, this scalar product allows 
for a generalization of the method of Ref. 38] in the case 
of inhomogeneous situations. 

If one considers solutions of the Schrodinger equation 
with a finite range potential as functions of r and of 
the energy E: ^[E,r), then a bound state of energy 
Eb and wave function b = '^{Eb,?) corresponds to 
a case where ^'(i?, r) is square integrable. Using the 
Schrodinger equation, one can show that: 



r<ro 



d^n{'<SJ*{E,r)drdE-^{E,r) 

Jr=ro 

-dr^*{E,r)dE^{E,r)) . (60) 



Using this identity permits to deduce that the normal- 
ization of a bound state is linked to the residue of the 
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scattering amplitude at E = Eb [s^- Using Ea. (|60|) . the 
norm of the state ^'(£',r) can be also written as: 

Jr>R 2^ Jr=B. 

. (61) 

The link between Eq. H61|l and the result given by the reg- 
ularized scalar product appears by noticing that Eq. H60|l 
is exactly the opposite of the r.h.s. of Ea. (|54l) in the 
limit k' ^ k and for = R. As a result, the norm of a 
state in the zero-range scheme obtained through Eq. H56() 
coincides with the formal limit i? — > in Ea. l|61(l . 

In the absence of an external potential, wave functions 
of two body bound states of energy Eb = — ?i^K^/2/i in 
channel (l) are of the form (for r > R): 

q,B{r,n)^M<l>imr'(^ldr^ (^^^^y^^ , (62) 

where (/);(r2) is a normalized angular function. The ex- 
pression for the bound state energy is found by solving 
1///(*k) = 0. For a partial scattering amplitude given by 
Eq. (|3J| , one finds in the resonant regime i. e. for large and 
positive values of , a bound state of vanishing energy 
with ~ 1/wiai (obviously in the case of the s-wave 
broad resonance ao = and k = I/wq). This is the state 
which plays a crucial role in the BCS-BEC crossover. In 
the zero-range approach, Ea. l|62|) is formally extended in 
the domain < r < R and using the regularized scalar 
product, one finds that for (\E'b|^'b)o — 1: 

This result is very interesting as it coincides with the one 
given by the residue of // at the energy E = Eb i38j. 

To conclude this section, we show that the scalar prod- 
uct gives a constraint on the possible values of the 
"effective range" parameter ai compatible with a consis- 
tent description of the low-energy bound state and hence 
of the resonant regime. Indeed, the probability of find- 
ing the molecular state outside the potential range is less 
than one: /^^^d^r I^Psp < 1. Then, using Eq.ll^SI) and 
keeping the dominant contribution as i? ^ 0, one finds 
that the parameter a; verifies in the resonant regime for 
I > 0: 

> (2;-3)!!(2Z- 1)!! . (64) 

In the case of a narrow s-wave resonance, in the regime 
of intermediate detuning {oq 3> wq) and wq > i?, then 
~ 1/woao and the condition uq ^ R follows directly 
from kR ^ 1 . By the way, Eq. (|64|l implies also that 
cxi > 0. As an example, the experimental fit performed in 
[29] for a particular p-wave resonance is compatible with 
Ea. (|64|l . Let us note that the ideal case for a description 
of the bound state within the zero-range approach corre- 
sponds to the situation where the two terms in Ea. 164(1 
are almost equal. In this regime the probability of find- 
ing the dimer in the region r > i?, i.e. in the region 
described by the zero-range model, is maximum. 



IV. CONCLUSIONS 

In this paper, we have developed a general zero-range 
one channel model for resonant scattering in arbitrary 
partial waves. In this approach, the wave function is a 
solution of the free Schrodinger equation and the inter- 
action is replaced by energy independent boundary con- 
ditions on the wave function as the distance between the 
two interacting particles goes formally to zero. This way, 
we generalize the Bethe-Peierls approach which has been 
widely used in the case of broad s-wave resonances. We 
have shown, how the boundary conditions can be imple- 
mented in the Schrodinger equation through zero-range 
pseudopotentials. In a given partial wave channel, the 
interaction can be described either in terms of a bound- 
ary condition on the wave function or with a family of 
pseudopotentials generated by an extra degree of free- 
dom (the parameter A). By construction, the Schrodinger 
equation is invariant under a change of the parameter A. 
This transformation appears then as a general symme- 
try of zero-range potential approaches. The pseudopo- 
tentials derived in this paper are the generalization of 
A-potential already obtained in the s-wave and p-wave 
channels . This class of pseudopotentials can serve 

as source terms in the Schrodinger equation permits us 
to search solutions with the Green's functions method. 
While exact treatments are invariant under a change of 
the A parameter, this property is in general no longer 
true in approximate schemes. This gives in turn a way 
to improve the approximations made (this idea can be 
applied in mean field treatments and has been already 
used in the Hartree-Fock-Bogoliubov formalism for the 
Bose gas 0,0)- For example, we have shown that with 
a specific choice of A, the first order Born approximation 
gives the exact result for two-body scattering states at a 
finite colliding energy. 

Using the fact that the scattering states are not mu- 
tually orthogonal in the zero-range scheme (except for 
the Fermi pseudopotential), we have shown how the no- 
tion of a regularized scalar product emerges naturally. 
Normalization of bound states by means of this metrics 
gives in free space the same result as the method based 
on the analyticy of the scattering amplitude. As it ap- 
plies directly in configuration space, this new tool gives 
us a simple way to perform the normalization in inhomo- 
geneous situations. 

The formalism presented in this paper can be used 
to obtain the scattering amplitude of particles confined 
in linear or planar atomic waveguides, thus allowing for 
studies of ultracold gases trapped in low dimensions. The 
few-body problem is another direction of research where 
this approach should be fruitful. 

Note added: Recently we learned of recent related 
works by Idziaszek and Calarco ^3 and also Derevianko 
psf where the form of the pseudopotential in Eq. (|53l) is 
also found. The main difference of the present approach 
with respect to Ref. |^ (or also Ref. 37]) is that the 
pseudopotential in Eq.J^HJ or equivalently the general- 
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ized Bethe-Peierls conditions in Ea. H16|l have been de- 
rived in an energy independent form. 



where (f2 = (0, 0), ' = {9',(j)')) are the angles defining 
the spherical coordinates for the two vectors (r, r '). The 
projection of a wave function ^ reads: 
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APPENDIX: USEFUL PROPERTIES OF THE 
SYMMETRIC TRACE FREE TENSORS Si 

This appendix collects some useful properties of the 
symmetric trace free tensors Si defined in Ea. H18(l lil . 
149] . Let us consider the two vectors f = x°'ea and f' = 
x'^Ca- By definition: 



(21-1) 



\r-r'\ l\ 

1=0 



-<Sllal3... 



(65) 



Introducing the unit vectors n — r/r = n°'ea, n' = 
r' jr' — n'°'ea and the symmetric trace free tensor S[ 
associated with f' we have: 



Si[ap...]n n 
so that Ea. (|65|l can be rewritten 
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1=0 



II 



(66) 



(67) 



This expression can be identified with the standard ex- 
pansion (for r ' < r) [H^ : 



1 1 / r 



r — r ' r ^ — ' \ r 



(68) 



where Pi{x) is the Legendre polynomial of degree I and 
we obtain: 

Viin.n')^^^^Si[^,,,,^S't^---^ . (69) 
For example, taking n' = shows that: 



Si{zzz...z) {21— ly}^^^ 



cost 



(70) 



with 9 = (f, e^). Ea. (|^ can be used to have an expres- 
sion of the projection operator 11; over the partial wave 
I in terms of the tensor Si. For that purpose, we recall 
the addition theorem 15(1: 



Vi{n.n') 



An 



21- 



_ r'™(i7)y'™*(f]') 

m— — l 



(71) 



(f|n;|*) = (2^ + 1) /^*(r,^2')^;("-"') • (72) 



together with Eq. (|69|l and the definition of the multipoles 
of the wave function Ea. H20|) . we obtain Ea. (|21|) . 

In this paper, we also use an interesting property of 
the symmetric traceless tensors found in Ref. [i^. Let 
us consider two symmetric and traceless tensors A and 
B of order I, then 



II 



{21- 



—A^^p,,,^B^-P-\ . (73) 



As an example, taking A = Si and B = S'{ allows one to 
recover from Eas. H69l73|) the known relation: 



Vi{n.n") = (2^ + 1) [f^Vi{n.n')Vi{n'.n"). (74) 



We close this appendix with the link between the sin- 
gularity of the Ith partial wave of a wave function with 
the derivatives of the delta distribution. To this end, we 
use the so-called relation: 



(75) 



together with Ea. l|?)7jl and the formal expansion: 



Sir-f') = J2'-^{dlp,,,^6){f) 



(76) 



where (9|^^ ]^){^) is the Ith partial derivative of the 
delta distribution (see Ea. l|19|l ). Finally, we obtain: 



Air t}^'lj dl, ,<5)(r)n'"n'^... . (77) 



{21-1 



^{d\^,,„]S){r)nV... . 



This expression justifies the introduction of the potential 

Vi^''' in Eq.(|2SI). Using Eq.JHHIl, one can also write Ea. l(77|) 
in terms of Legendre Polynomials: 



Viicost 



= 47r 



(-1) 



(78) 
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